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Q\ ' Abstract 

, We study the process e^e+ — > X1X2 ^i^h the subsequent decay X2 ~^ Xi^^^^ taking into account the complete 

O 1' spin correlations. We give the analytical formulae for the differential cross section and present numerical results for 

I , the lepton angular distribution and the distribution of the opening angle between the outgoing leptons for the LEP2 
energy of ^/s ~ 193 GeV and for = 500 GeV. We examine three representative mixing scenarios in the MSSM 

(-H and also study the influence of the common scalar mass parameter mo on the shape of the distributions. For the 

• • i lepton angular distribution the effect of the spin correlations amounts to up to 20% for lower energies. Wc find that 

, ^ the opening angle distribution is suitable for distinguishing between higgsino-like and gaugino-like neutralinos. The 

^ , shape of the lepton angular distribution is very sensitive to the mixing in the gaugino sector and to the value of mo. 

yu, For higher energies it is also suitable for distinguishing between higgsino-like and gaugino-like neutralino. 
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I. INTRODUCTION 



Although the Standard Model (SM) is extraordinarily successful for describing the electroweak phenomena, it has 
several theoretical shortcomings. The most severe of these deficiencies, the hierarchy problem, can satisfactorily be 
solved by the concept of supersymmetry (SUSY) broken at the TeV scale. This concept received a fresh impetus by 
the result that in a SUSY-GUT model the measured coupling constants evolved to high energies meet at a single 
point. 

The most economical candidate for a realistic SUSY model with minimal gauge group SU{3) x SU{2) x U(l) and 
with minimal content of particles is the Minimal Supersymmetric Extension of the Standard Model (MSSM). In this 
paper we consider its simplest version with the conserved quantum number R parity. Its implications include that 
SUSY particles can only be produced in pairs and the lightest supersymmetric particle (LSP) is stable and escapes 
detection. As usual we assume that this particle is the lightest neutralino Xi- 

Among the new particles the charginos (the supersymmetric partners of the charged gauge and Higgs bosons) and 
the neutralinos (the partners of the neutral gauge and Higgs bosons) are of particular interest. As they are expected 
to be lighter than the gluino and in most scenarios lighter than the squarks and sleptons , [§) , the next-to- lightest 
neutralino X2 ^^'^ the lightest chargino could be first observed in future experiments at e^e^-colliders. In particular 
the production of X1X2 pahs allows to study a wide region of the SUSY parameter space. Although in general chargino 
production is favoured by larger cross sections, in certain regions of the parameter space sizeable cross sections for 
the neutralino process can be expected ^ . Moreover it might be possible to discover SUSY by neutralino production 
if charginos are not accessible. 

If new particles are discovered which are possible neutralino candidates, for a clear identification the complete 
investigation of their decay characteristics is indispensable. Neutralino decay widths and branching ratios have 
been thoroughly studied Q|. Angular distributions and angular correlations of the decay products can give valuable 
additional information on their composition from photino, zino and higgsino components. 

Another interesting question is to see if angular distributions of decay products allow separation of neutralino 
production from chargino production. 

Moreover, from decay angular distributions with complete spin correlations of the decaying particle one can infer 
the spin of the new particles. 

Finally the identification of neutralinos would be completed by ascertaining their Majorana character. In Ref. 
it is demonstrated that this is possible by means of the energy distributions of the decay leptons if the neutralinos 
are produced in collisions of polarized e^e~ beams. The angular distributions of the decay products might, however, 
offer the possibility to prove the Majorana character if polarized beams are not available. Furthermore the angular 
distributions of the final leptons are suitable observables for studying CP-violation in supersymmetric models 0,0. 

The above mentioned reasons motivate a study of angular distributions in associated production of neutralinos and 
the subsequent decay of the next-to-lightest neutralino. Since angular distributions depend on the polarization of the 
parent particles one has to take into account spin correlations between production and decay. 

In general quantum mechanical interference effects between various polarization states of the parent particles pre- 
clude simple factorization of the differential cross section into a production factor and a decay factor |^ unless the 
production amplitude is dominated by a single spin component However, in Ref. it is shown that the factoriza- 
tion property holds for particles with spin provided that suitable Lorentz invariant variables are used. But this is not 
the case for angular distributions of the neutralino (or chargino) decay products in the laboratory frame. For energy 
distributions of the final particles in the laboratory frame the spin correlations are nevertheless usually ignored. 

Heavy fermion production with subsequent decay was considered in Ref. [|lO|. Recently a Monte-Carlo generator 
for chargino production and decay including spin correlations was developed in Ref. . 

In this paper we calculate the cross section of the process e^e^ — s- X1X2 ^^'^ subsequent direct leptonic decay, 
X2 ~^ We give the analytical formulae of the differential cross section with complete spin correlations of the 

decaying neutralino. We study numerically the influence of these spin correlations on energy and angular distributions. 
The energy distributions are independent of spin correlations. However, for the lepton angular distributions for lower 
energies the effect of the spin correlations amounts to up to 20%. The shape of the lepton angular distribution is 
very sensitive to the mixing character in the gaugino sector and to the value of ttiq. For higher energies it is suitable 
for distinguishing between higgsino- and gaugino-like neutralino. We also analyze the opening angle distributions 
and show that at lower energies they are suitable for distinguishing between higgsino- and gaugino-like neutralinos. 
This will allow one to constrain the parameter space of the MSSM. We also consider the influence of the scalar mass 
parameter mo on the shape of distributions. 

In Sec. 2 the general formalism is shown. In Sec. 2.1 the Lagrangian, couplings and Feynman diagrams are given, 
in Sec. 2.2 the spin-density formalism is presented, and in Sec. 2.3 the kinematics is given. In Sec. 3 we present the 
formulae for the differential cross section with complete spin correlations. Numerical results for the LEP2 energy 
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y/s = 193 GeV and for -^s = 500 GeV and a discussion are presented in Sec. 4. 



II. GENERAL FORMALISM 
A. The Feynman diagrams 

In this section we show the Feynman diagrams and give the Lagrangian for the production process, e~ {ki)e^ {k2) — > 
Xi{<li)X2{l2) , and for the direct leptonic decay, X2(92) x\{Pi)^'^ {P2)i~ ijps.)- The arguments ki,k2,qi,q2 and 
Pi)P2,P3 denote the momenta of the incoming electron, positron, the produced neutrahnos X11X2 and the outgoing 
neutrahno Xi and leptons i'^ , £~ from the X2 decay. Both the production and the decay process contain contributions 
from Z'^ exchange in the direct channel (s-channel) and from cl and en exchange in the crossed channels (t-, it- 
channel). We introduce the kinematic variables: 



(fci- 

(fc2 
(fcl 



k2r, 
92)^ 



{P2 
(92 - 
= (92 



-P3) 
P2)- 

-Ps) 



(1) 

(2) 
(3) 



Channels referring to the decay are marked by a dash. The Feynman diagrams are shown in FIG.|l|. 
From the interaction Lagrangian of the MSSM ( in our notation and conventions we follow closely 



£z<> 



zf>e+e 



cos 9w 



Z,Jj^[LePL + RePR]e + h.c. 



yw 



J^uio = gfu^PRX'lh + gf^UPLxT^B^ + h.c, J = 1, • • • , 4 



(4) 

(5) 
(6) 



we obtain the couplings: 
• Li ^ T^i - a sin^ 9w, Ri = -ee sin^ 9w 



-V2 



cos t)\Y 

^ — —V^ee sin6'H' 



3i ■ 



ei sin^ 9w)Ni2 



€£ sin 9wNii 



tan 9wN*^ 



^13 



■ \ (n^N*^ - N,^N*^ cos 2/3 - i (n,^N*^ + N,^N*^"^ sin 2/3, 



In our case i = 1, j = 2 and we shall write = O^'^. Here Pl,r = \{X T 75)7 9 is the weak coupling constant 

{g — e/ sin^vi^, e > 0), and ei and T^i denote the charge and the third component of the weak isospin of the lepton 
I. Furthermore tan/? = vi^2 are the vacuum expectation values of the two neutral Higgs fields and Nij is the 

unitary 4x4 matrix which diagonalizes the neutral gaugino-higgsino mass matrix in the basis 7, Z, H^. Since we 
disregard in this paper CP violating phenomena the elements Nij of the diagonalized matrix and the couplings can 
be chosen real. Then some of the neutralino mass eigenvalues may be negative and we shall write them in the form 
rjim.i with rrii > 0,77^ — ±1 |13|. The respective amplitudes for the Feynman diagrams are taken from Ref. fl^. 



B. Spin density matrix for production and decay 

The differential cross section for the production e^(/ci)e+(A;2) ^ Xi(9i)X2(92) and the subsequent decay X2(92) ^ 
XiiPi)^'^ {P2)(r {p-i) is given by 

da = ^ \ ^ j r,"^Yr^'l'(2^)'^^'Hfci +k2-qi - Pi - P2 - P3)dlips{q,,p,,p2,P3), (7) 

where s is the cms-energy squared of the incoming e~ and and X{x,y,z) — x'^ + y'^ + z'^ ~~ 2xy — 2xz — 2yz is 
the kinematical triangle function. All lepton masses are neglected. In Eq. (7) Ae-,Ae+,Ai and Ai,A+, A- label the 
helicities of the electron, positron and xf of the production process and the hehcities of the Xi, of the decay 
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process, respectively. dlips{qi,pi,p2,p^) = (2^)3217" ' ' ' (27r)3'2p" Lorentz invariant phase space element of the 

four final particles. The helicities A^- , \^+ of the initial particles are averaged because in this paper we consider the 
case of unpolarized beams. Since the polarization of the outgoing particles will not be measured, the helicities of the 
final particles Aj and Ai, A-)-, A_ have to be summed over. Therefore we suppress these helicity indices in the following 
and only display the helicity A2 of x"- 

The total width T2 of the decaying neutralino X2 is small compared to its mass m2 . Therefore the amplitude T of 
the combined process is a coherent sum over all polarization states of the helicity amplitude P^^ for the production 
process times the helicity amplitude for the decay process and a pseudopropagator A2 = -^ri — -rT-- — p- of X2- 

T = A2P^^Dx, (8) 

The amplitude squared for the combined process 

|Tp = |A2|V^^"vf,,, (9) 



is thus composed from the unnormalized spin density matrix 
of the neutralino X2 ^-nd the decay matrix 

Py,x,^D,,Dl, (11) 

for the respective decay channel. As in FIG|| all helicity indices but that of the decaying neutralino are suppressed. 
Repeated indices are summed over. 

Squaring the total amplitude one obtains interference terms between various helicity amplitudes. These terms 
preclude factorization in a production factor jP'^'^p times a decay factor X]A2l^A2p as for the case of spinless 
particles. 

We use the general formalism to calculate the helicity amplitudes for production and decay of a particle with four- 
momentum p and mass m. Therefore we introduce three spacelike four-vectors s^, (a = 1,2,3) which together with 
p/m form an orthonormal set: 

(12) 

m 

s'^-s'' = -S"'' (13) 
sl-4 = -9,^ + ^- (14) 

A convenient choice for the explicit form of s° is in a coordinate system where the direction of the three-momentum 
is p (sin O cos $, sin sin<I>, cos 0) | ]T5| : 

s^^ = (0,-sin$,cos$,0) (15) 
s^M ^ (0, cos 9 cos $, cos 9 sin $, - sin 9) (16) 

s^>^^-{\p\,Ep). (17) 

TO 

Then in this frame of reference s(i)^(^) and s'^' describe transverse and longitudinal polarization of the particle. 

When computing the density matrices for production and decay, Eqs. (10), (11), we make use of the Bouchiat-A/Iichel 
formulae Jisf : 

u{p,X')u{p,X) = i[<5AA' +75 raly]{^ + m) (18) 



v{p,X')v{p,X) = i[<5AA'+75 /V^,J(^-m). (19) 

In the amplitude squared (Eq. (9)) the spin vectors s" linearly enter the matrices pp (Eq. (10)) and (Eq. (11)). 
They induce by Eq. (14) the above mentioned quantum mechanical correlations between production and decay. 
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C. Kinematics and Phase Space 



We split the phase space into the one for the production, e^(fci)e+(fc2) — ^ Xi(9i)X2('Z2), and the one for the three 
particle decay, X2(92) ~> Xi(pi)^^(P2)^ (ps)- Then we obtain the differential cross section in the e~e'^~cms by 
integrating over the effective mass squared S2 of the decaying neutralino X2 101 ■ 



da 



ds2 W 



8E^ J 2tt (s2 - "7,2)2 + r2m2 



{2'kY2Ei {2tj:Y2E2 



^ (o ^4x(4)/ , , ^ d^pi d^p2 d^p3 



with W — pp ^ Py \2 (compare Eq. (9) ). The energies of the produced X17X2 are Ei,E2 and those of x\,t^ and I 
from the decay are denoted by Ei , E^ and . Ei, is the beam energy, m2 is the mass of X2 and the mass of Xi is 
labeled by mi. Finally 6m- (0m+) is the angle between the electron beam and the outgoing negatively (positively) 
charged lepton t~ (^+) in the laboratory frame (FIG.||). 

If the total width of the decaying particle X2 is much smaller than its mass, r2 <C m2 , one can make the narrow width 

approximation -, ri- — — ^(5(s2 — rnV\ and one obtains for the differential cross section in the e"'"e~-cms, 

i.e. in the laboratory frame: 

1 1 \q2\ E+E^dE^sinOded^dn+dn. 

16(2^)7 \m2\r2 64£;3 [E2 - Igal cos 62+ -E^{1- cos e+_)] ' ^ ^' 

Q is the production angle of X2; ^ describes the rotation around the electron beam axis, and dfl+{dil^) is the solid 
angle of £+{£-). 

In consequence of momentum conservation the energy E^ is determined by the angles 02- (^2+) between i~ 
(£+) and X2j and the opening angle 6-| between the leptons: 

^ m2-m2-2j;_(ig2-|g2|cos02-) .221 

+ 2[£;2-|g2|cos6l2+-£;-(l-cose+_)]' ^ ' 

The possible region of the lepton energy E- depends on the lepton decay angle 92-- 



so that E+ = 0, see Eq. (22), at the upper bound. 



III. ANALYTICAL FORMULAE 

In this section we give the analytical expressions for the product W ~ Pp ^^PvAa (compare Eq. (9)) of the den- 
sity matrices for the production process, e~{ki)e~^{k2) Xi(9i)X2(92), and the direct leptonic decay, X2(92) 
XiiPi)^^ {P2)£~ {Ps) ■ Both the production and the decay process contain contributions from Z° exchange in the direct 
channel and from and exchange in the crossed channels (compare FIG|^). 

The product W in Eq. (9) is a sum of contributions from the different production and decay amplitudes and their 
interference terms, 

W=Y1 ^-b.cd, (24) 

where the first pair a, b of indices refers to the production process, and the second pair c, d refers to the decay process. 
The pairs (ab) and (cd) run through all combinations of the values Z, L4, L„, i?„, where Z and L,R denote the 
exchanged particles Z'^ and e^, en and £l, respectively. In the case of slepton exchange the indices t, u denote the 
channel. Thus WzLt,RtR^ is the contribution of the interference term to the production process from Z°-exchange 
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and e^-exchange in the t-channel and the interference term to the decay process from ^^j-exchange in the t-channel 
and £i{-exchange in the u-channel (FlCjl]). 

There are altogether 121 contributions Wab,cd which can be classified into 16 groups. We give one representative of 
each group and list the indices of the possible combinations. The other terms Wat.cd of this group are then obtained 
by substituting the propagator combination according to the pair of indices (see Eqs. (65)-(70)) and by substituting 
momenta and couplings and sign factors as given in TABLE I below. The substitutions of the momenta ki,k2,P2, Ps 
and of the couplings 0^,0^,L,R also have to be performed in Eqs. (41)~(64). In the representative term the sign 
factors ^,v,ijj,T,v have the value +1. 

For illustration we consider as an example WLtL^,LtLti which is the representative of Group 12. In order to 
get the contribution of Wr^r^^l^l^ one has to change the propagator combination AL(t)AL(M)AL(i)Ai(f) 
A7j(t)Afl(u)A2^(u)Ai(M) and one has to substitute cr in the production process and one has to substitute 

t — *■ u in the decay process. When substituting cl — >■ eR, one has to change the sign factors i^, v and the couplings 
according to the first line of TABLE I. Moreover, when substituting t u one has to change the sign factors t, v, the 
couplings and the momenta p2 Ps according to the fourth line of TABLE I. 

Group 1(1 term): 

Wzz.zz - 32|Az(s)nAz(s)p|i?l • PI - ■ 5I2 + DI2 ■ SI3} (25) 



Group 2 (4 terms): 
Production: {LtZ), {RtZ), (L^Z), (R^Z) Decay; {ZZ) 



Wl,z,zz = 16AL(i)| Az(s)p{ - fiRe{Azis)) [Dl ■ P2 + v{D23 ■ 522 - D22 ■ S23)] 

+ i^vIm{Az{s))[D23 ■ S62 - D22 ■ 563]} (26) 



Group 3 (4 terms): 
Production: {LtLt),{RtRt),{LuLu),{RuRu) Decay: {ZZ) 



Wl,l,,zz = 8Ai(i)|Az(s)p{(A:igi)(fc2g2)/?l - v{D23 ■ S?,2 - D22 ■ SSj)} (27) 



Group 4 (2 terms): 
Production: {LfLu), (i?ti?„) Decay: (ZZ) 



Wl,l^,zz = -8AL(i)AL(7i)|Az(s)p{Dl • P4 + v7^imi{D2i ■ SA2 - D22 ■ SA^)^ (28) 



Group 5 (4 terms): 

Production: {ZZ) Decay: {LtZ), {RtZ), {LuZ), {R^Z) 

WzzMZ = 16|Az(5)pAi(t)r{i?e(A,(s)){Pl • D3 - v[DA ■ SI2 - D5 ■ Sh]} 

- Im{Az{s)){Rl - Ll)0^rjimi[S72 ■ P3i - S7i ■ PSa]} (29) 



Group 6 (16 terms): 

Production: {LtZ), {RtZ), {L^Z), {R.,Z) Decay: {LtZ), {RtZ), {L^Z), {R.,Z) 

WL,ZX,Z = 8AL{t)AL{t)T 

■ I - iiRe{Az{s))Re{Az{s))[P2 ■ D3 + v{D4: ■ 522 - -05 • 523)] 

+LuvIm{Az{s))Re{Az{s))[D4, ■ S62 - -D5 • 563] 
~iyRe{Az{s))Im{Az{s))Leriimi 

■ [2772m20^(A:igi)572 + mmiO"^ {hq2)S72 - ryi"^iO^(fc292)57i] 
+Im{Az{s))Im{Az{s))L,0'^L,0'^ml[{p2q2)S'i3 - (^392)^42 + mlS5]} (30) 
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Group 7 (16 term,s): 

Production: (LtLt), (RtRt), (LuLu), (RuRu) Decay: (LtZ), (RtZ), (L^Z), (RuZ) 

Wl,l„l,z = 4Ai(f)Ai(t)r{i?e(Az(s))[(A:igi)(A;2g2)£'3 - v{D4 ■ 3^2 - Db ■ SSs)] 

+ a;7m(Az(s))r?imir/2m2(fcigi)S'72| (31) 

Group 8 (8 terms): 

Production: (L^L^), Decay: {LtZ),{RtZ),{L^Z),{R^Z) 

Wl,l^,l,z = -4Ai(i)Ai(u)Ai(t) • T[Re{Az{s))[D3 ■ PA + vr]imi{D4: ■ S42 - D5 ■ S^a)] 

+ ojIm{Az{s))ml[{hq2)S72 - {k2q2)S7i]} (32) 

Group 9 (4 terms): 

Production: (ZZ) Decay: {LtLt), {RtRt), (LuLu), {RuRu) 

WzzMU = 8|Az(s)pAi(i)(pip3) • {(P2?2)P1 + u??2m25l2} (33) 

Group 10 (16 terms): 

Production: (L^Z), {RtZ), (L^Z), (R^Z) Decay: (LtLt), (RtRt), (i„i„), (RuRu) 

Wl,z,l,l, = 4Az,(t)Ai(i)(pip3) • { - iJ.Re{Az{s))[{p2q2)P2 - vr]2m2S22] - uvIm{Az{s))'n2m2S62} (34) 

Group 11 (16 terms): 

Production: {LtLt), {RtRt), {LuLu), {RuRu) Decay: {LtLt), {RtRt), {LuLu), {RuRu) 

Wl,l„l,l. = 2Al{t)Al{i)[{pips){p2q2){krqi){k2q2) + vr}2m2{pWi)S2,2] (35) 

Group 12 (8 terms): 

Production: {LtLu), {RtRu) Decay: {L,Li). (RiRt), {L^L^), {R^Ru) 

Wl,l^,l,l, = -2Ai(t)Ai(w)A|(t){(piP3)(f^2g2)i'4 - vr]imir]2m2{piP3)S42} (36) 

Group 13 (2 terms): 
Production: {ZZ) Decay: (LfLu), {RtRu) 

WzzMLu = 8|Az(s)pAt(t)Az,(u){r;iTOir;2m2(j?2P3)Pl + vr)imi[{pzq2)Sl2 - (j32g2)5'l3]} (37) 

Group 14- (8 term,s): 
Production: {LtZ),{RtZ),{LuZ),{RuZ) Decay: {LtLu),{RtRu) 

Wl.zml^ = '^AL{t)AL{t)AL{u) ■ { - iJLRe{Az{s))[r^imir]2m2{p2P3)P2 + vr^imi{-{paq2)S22 + {p2q2)S23)] 

+ vvIm{Az{s))iqimi[-{piq2)S&2 + (^292)^63]} (38) 

Group 15 (8 terms): 
Production: {LtLt), {RtRt), {LuLu), {RuRu) Decay: {LtL^,), {RtRu) 

WLtLt,LtL„ = 2Al{t)AL{F)AL{u)r]imir]2m2{kiqi) ■ ^{p2P3){k2q2) + v[-{p3q2){k2P2) + {p2q2){k2P3)]^ (39) 
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Group 16 (4 terms): 
Production: (LtLu), {RtRu) Decay: (LtLu), (RtRu) 

Wl,l„,l,l^ = -2AL{t)AL{u)AL{t}AL{u)r]nni ■ {v2m2{p2P3)P'i + vr]imi[-{p3q2)S42 + {p2q2)S43]^ (40) 

For the sake of a clear presentation of our analytical results we have introduced three groups of abbreviations. The 
first group refers to the production process: 

PI = [Rl + LDiO"^" [ihcn){k2q2) + (fcig2)(fc2<Zi)] + 0'^0%mi772m2(fciA;2)} (41) 

P2 = 2L,0^{kiqi)ik2q2) + mmimm2LeO^{kik2) (42) 

P3i = mmiO^{kiq2) + r]2m20^{kiqi) (43) 

P32 = VimiO^{k2q2) + mm20^{k2qi) (44) 

P4 = r]imirj2m2{kik2) (45) 

and the second group refers to the decay process: 

Dl = {Rl + LI)0^ {{piPz){P2q2) + (piP2)(P392) + r]i'mir]2m2{p2P3)} (46) 

£122 = {Rl - i^)0^'{r?imi(p2(72) +r?2m2(piP2)} (47) 

D23 = {Rl - Ll)0''\r]imi{p3q2) + V2m2{piP3)} (48) 

D3 = ie[20^(piP3)(P292) - 0%imir]2m2{P2P3)] (49) 

DA = L„[dSimi{p3q2) - 20^?72m2(piP3)] (50) 

D5 = LeO%imi{p2q2). (51) 

The third class is related to the spin correlations between production and decay and connects the momenta of both 

subprocesses: 



512 = {Rl - i^){?]imiO«0^[(fc2P2)(fci<Z2) - {kiP2){k2q2)] 

+mm20-'[ik2,M- {kW2) + ihS^liP^) - {k,qM- {k2P2) + (^M^hM)]} (52) 

513 = {Rl - Ll){r,imiO''0''[{k2P3){kiq2) - (fciP3)(fc2g2)] 

+,2^20-^ [{k2q,) ( - {k,ps) + ihSllM) _ ( _ ^k2Ps) + (^^^2^) } (53) 

(^2g2)(P2«2), 



^2 



522 = -2j?2m2LeO^(A;igi)[-(fc2P2) + 

+mmiLeO^[{k2P2){kiq2) - {kiP2){k2q2)] (54) 

523 = -2r^2m2L^O^{k,q,)[-{k2P3) + i^Sl}^] 
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+VimiLeO^[{k2p:i){kiq2) - {kip3){k2q2)] (55) 

CQ fl, W fl \ , (^292)(P2g2)i 

S32^V2m2{kiqi)[-{k2P2)^ 2 J (56) 

533 = V2m2{km)[-{k2Ps) + i^H^^] (57) 



m2 

542 = {k2P2){kiq2) - {kip2){k2q2) (58) 

543 = {k2P3){kiq2) - {klP3){k2q2) (59) 

55 = {kip3){k2P2) - {kiP2){k2P3) (60) 

562 = Le0^riimi[k2kip2q2] (61) 

563 = Le0^riimi[k2kiP3q2] (62) 

571 = IeO^[(Z2fclP3P2] (63) 

572 = LeO'^[q2k2P3P2] (64) 
with [abcd\ = e^vp„a'^b"cP(F . 
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These terms, i.e. Eqs. (52)-(64), would be missing if we had assumed factorization of the differential cross section 
in production and decay. 

For a better transparency all couplings originating from the decay process are marked by a dash. The indices are 
used in order to emphasize the symmetry between and e+ in the initial state and and in the final state, 
respectively. 

We have introduced the following products of propagators and coupling constants: 

^^(^) = 2 ! ■ r (65) 

-/^^^ , ^n{t)^j^-fg*-fg (66) 

AUu) = ^^-ft-f[,* , An{u) = ■ ff, ■ f,f (67) 

u — rat u — mi 

ei, en 

cos^ dm s — + irnzL z 
^^(^") = f:^-/^^2*-#i , Anii) = ■ fg* ■ fj, (69) 

El 

AL{u) = ^^^-fg-fn , An{u)^^^^-fg-f,f. (70) 

u — mj u — mi 



and 



IV. NUMERICAL RESULTS AND DISCUSSION 



A. Scenarios 

Neutralinos are linear superpositions of the photino 7, the Zino Z and the two higgsinos and H^. The 7 and 
Z components only couple to the selectrons whereas the higgsino components couple to the Z^. The composition 
and the masses of the neutralino states depend on the three SUSY mass parameters M, M' (sometimes also called 
M2 and Ml) and /i, whose values follow from the specific SUSY breaking mechanism, and on the ratio tan/3 = vijvx 
of the vacuum expectation values of the Higgs fields. In order to reduce the number of parameters we shall assume 
M' — %_Mtw?Bw as suggested by grand unification The gaugino mass parameter M is related to the gluino 
mass by M R:! O.Srrig and the gluino mass is roughly given by m-g « 2.4m]^/2 P^ ]? where TO1/2 is the common 
gaugino mass at Mcut(m ~ 0.72mi/2)- 

The masses of the sleptons are related to the SUSY parameters M and tan f3 and to the common scalar mass 
parameter mo at the unification point [ pTf : 

m| = ml + 0.79A/2 + ml cos 2/3(-0.5 + sin^ Ow) (71) 

= Too + 0.23M2 - ml cos 2/3 sin^ Ow (72) 

In order to illustrate the influence of the neutralino mixing and of the scalar mass mo we shall consider three 
representative scenarios which differ significantly in the nature of the two lowest mass eigenstates Xi an.^ X2- The 
selectron masses are calculated for two values of the scalar mass mo, mo = 90 GeV and mo = 200 GeV. For the 
parameters of the Standard Model (SM) we take mz = 91.19 GeV, Tz = 2.49 GeV, sin^ 6*^^ = 0.23 @ and 
a — 1/128. We choose tan/3 = 2. The parameters of our scenarios and the mass eigenvalues of the two lightest 
neutralinos, the light chargino and the selectrons are given in TABLE II. The width of the X2 has been computed 
with the program of [|§ (see TABLE II). 

Notice that also in scenarios with /i < the branching ratio for the radiative decay is less than 0.5% in the examined 
region of parameter space ||2|] . 

The 7 and Z are mixtures of the B and the W3 gauginos, 7 — cos 9wB + sinOwWs, Z — — sin 9wB + cosOwWs. 
Therefore in TABLE III the components of the neutralino states are given in the basis (7, Z, H^, H^) and in the basis 

(B,H^3,^°,i?,°). 

In scenario A Xi has a dominating photino component and X2 has a dominating zino component, whereas in scenario 
B both neutralinos are nearly equal photino-zino mixtures. In the two scenarios A and B the Xi is almost a pure 
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B-ino and the neutralino X2 nearly a pure Wa-ino. In scenario C both neutrahno states are dominated by strong 
higgsino components (TABLE III). 

For the phase space integration (chosen relative accuracy 10~^) we used the Monte-Carlo routine Vegas. The 
evaluation was made for cms-energies of y/s ~ 193 GcV and — 500 GeV. 

It can be derived from that the total cross section factorizes, i.e. it is independent of the spin correlations. This 
fact has been used in our numerical calculations as a check for the phase space integration and the total cross sections 
are given in TABLE IV. 

B. Lepton angular distributions 

In this section we give numerical results for the angular distributions of £~ with respect to the electron beam axis 
computed with complete spin correlations according to Eqs. (52)-(64). The angular distributions of £^ are obtained 
by substituting cosOm- — cos6a/_|_. 

In order to demonstrate the significance of the spin correlations we compare our results with those obtained from 
the assumption of factorization of the differential cross section into production and decay. As can be seen from 
FIGS.MHld for aU mixing scenarios and both values of the scalar mass mo the contribution of the spin correlations 
has the biggest effect in the forward and in the backward direction and vanishes in the direction perpendicular to the 
beam axis (cosGm- = 0). The contribution of spin correlations in the forward direction has opposite sign of that in 
the backward direction. Their magnitude decreases with increasing energy. 

Especially for lower energies the spin effect is sizeable in scenario A with a photino-like Xi and a zino-like X2 for 
both values of mo- For y/s — 193 GcV (^/s = 500 GeV) its magnitude amounts to about 20% (6%) in the forward 
and backward direction for both values of mo, mo = 90 GeV and mo = 200 GeV, FIG.|| (FIG.||). In both cases the 
contribution of the spin correlation is negative in the backward direction and positive in the forward direction. 

A comparison with the results for the gaugino-like scenario B (FIGJ| and FIG.^ shows how sensitively the spin 
correlation effect depends on the gaugino- higgsino mixing and on the value of mo- It is noteworthy that although in 
both scenarios A and B the x? is B-ino-like and X2 is Wa-ino-like with, however, different phases of the Wa-ino and 
B-ino admixture. For -^s = 193 GeV and mo = 90 GeV (FIG.^) the magnitude of the spin effect is only 2% in the 
forward and in the backward direction whereas for mo = 200 GeV (FIG.|^) it amounts to 10%. For higher energy 
yG = 500 GeV it is negligible (FIG.|). 

It is remarkable that in the case of gaugino-like neutralinos the scalar mass mo crucially determines the shape of 
the angular distributions. This is most obvious for scenario B and ^/s = 193 GeV. For the smaller value mo = 90 GeV 
the angular distribution has a maximum nearly perpendicular to the beam direction and is almost FB symmetric, 
ApB = —0.8% (compare TABLE V). The contribution of the spin correlations is positive in the backward direction 
and negative in the forward direction (FIG.^. Increasing the value of mo completely changes the shape of the angular 
distribution and for mo ~ 200 GeV it has a minimum in the backward hemisphere and the forward direction is 
favoured, ApB — -f5.9% (FIG.|^). Now the contribution of the spin correlations is negative in the backward direction 
and positive in the forward direction. 

In the higgsino-like scenario C both production and decay are dominated by Z'^-exchange. Therefore the dependence 
on mo is considerably smaller and we give only numerical results for mo = 90 GeV. Here the contribution of spin 
correlations is negligible and the angular distribution is practically flat for y/s = 193 GeV (FIG.^) and FB-symmetric, 
ApB = 0.07%, and for y/s — 500 GeV with a minimum perpendicular to the beam direction and ApB — 0.11% 
(FIGjlol). 

Notice however that a negligible FB-asymmetry is not an unequivocal signature for higgsino-like neutralinos. As 
can be seen from FIG.^ for scenario B with mo = 90 GcV the FB-asymmetry may also be small for gaugino-like 
neutralinos (see TABLE V). From FIGs.|, | and |l^ we conclude that for higher energies far enough from threshold 
the lepton angular distribution is suitable for distinguishing between gaugino-like and higgsino-like neutralinos. For 
y/s = 500 GeV the angular distribution is practically FB-symmetric with however a maximum for gauginos but a 
minimum for higgsinos perpendicular to the beam direction. 

C. The lepton opening angle distribution 

In contrast to the lepton angular distribution, the distribution of the opening angle between both leptons factorizes. 
Due to the Majorana character of the decaying neutralino the spin correlation terms are just cancelled by this partial 
phase space integration. 
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It is also noteworthy that in contrast to the lepton angular distributions, the distributions of the lepton opening 
angle are similar for both gaugino-like scenarios A and B (FIGsjl| and Especially for = 193 GeV they differ, 
however, distinctively from that for higgsino-like neutralinos in scenario C (FIG.p^). 

In the case of gaugino-like neutralinos the distributions for y/s = 193 GeV are rather flat (FIGs.|ll] and [L2h. 
Changing the scalar mass from mo = 90 GeV to toq = 200 GeV results in a reduction of the cross sections by a factor 
of approx. 3 in scenario A and by a factor of approx. 4 in scenario B. The shape of the distribution, however, remains 
essentially unchanged. 

One should note that for the shape of the opening angle distributions the influence of varying the value of mo is 
much smaller than for the lepton angular distributions. 

For higgsino-like neutralinos (FIG.p^), the shape of the opening angle distribution for ^/s — 193 GeV is completely 
different from those of gaugino-like neutralinos. Here the lepton pairs are preferably emitted with small angles between 
them, approximately 66% of them with an opening angle between and 7r/2. 

With increasing energy the opening angle distribution shrinks more and more and for ^/s = 500 GeV it displays for 

both higgsino- and gaugino-like neutralinos a rather narrow peak at or near cosB^ = 1. Obviously the distribution 

of the opening angle between the leptons is suitable to distinguish between gaugino- and higgsino-like neutralinos at 
lower energies. 

D. Energy Distributions 

Just as the opening angle distribution the energy distribution of the outgoing lepton factorizes. As a consequence 
of CP invariance (Sec. 2.1) and the Majorana character of the neutralinos the energy spectra of both leptons, (.~ and 
, are identical 

In FIG.|lJwe give the energy distributions for scenario B for mo = 90 GeV and mo = 200 GeV and cms-energies of 
y/s — 193 GeV and ^/s = 500 GeV. For all scenarios the position of the maximum is independent of the actual value 
of mo. 

E. Summary 

In this paper we have calculated the analytical expression for the differential cross section of the associated pro- 
duction of neutralinos, e~ -I- e+ X? + X2j subsequent direct leptonic decay, X2 ~^ Xi + + ^ i with 
complete spin correlations between production and decay. The angular and the energy distribution of the outgoing 
lepton as well as the distribution of the opening angle between both leptons have been computed for cms-energies of 
y/s = 193 GeV and yfs = 500 GeV. These distributions have been examined with regard to their dependence on spin 
correlations, on the neutralino mixing character and on the scalar mass parameter mo- 

The quantum mechanical interference terms between the various polarization states of the decaying neutralino give 
rise to a strong effect in the lepton angular distribution with respect to the beam axis, whereas the opening angle 
distribution and the energy distribution are independent from these spin correlations. 

For energies not too far above the threshold (y^ = 193 GeV) the opening angle distribution turns out to be 
suitable for distinguishing between higgsino-like and gaugino-like neutralinos. However, it is rather insensible to 
variable mixing in the gaugino sector. Here the shape of the opening angle distribution only slightly depends on the 
scalar mass mo- 

The lepton angular distribution, on the other hand, is for lower energies not only very sensitive to the mixing in 
the gaugino sector but also to the actual value of mo. 

For energies far above threshold the shape of the lepton angular distribution is rather insensible to the mixing in the 
gaugino sector and to the value of mo- It is, however, extremely different for gaugino- and higgsino-like neutralinos 
and suitable for distinguishing between them. 

If the neutralinos are gaugino-like, the effect of spin correlations in the angular distributions can be large amounting 
to about 20%. For higgsino-like neutralinos, on the contrary, the contribution of spin correlations is practically 
negligible. 

The energy distributions and the distributions of the opening angle, finally, are independent from the spin correla- 
tions. Apart from the magnitude of the cross sections they are rather insensitive to the actual value of mo. 

The clear structure of the analytical formulae presented here allows to include hadronic decays and to extend the 
investigations to cascade decays and to production and decay of, for instance, X2X2 pairs. Concerning the hadronic 
decays, X2 ~^ XiQQi where the outgoing quarks develop two jets, we expect a similar shape of the opening angle 
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distribution between quark and antiquark. Thus, the two jets would be better separated for gaugino-like neutralinos 
than for higgsino-hke neutrahnos which prefer small opening angle. 

These investigations as well as the inclusion of beam polarization and the results for the chargino process angular 
distributions will be discussed in a forthcoming paper with regard to the determination of SUSY parameters. 



V. ACKNOWLEDGEMENT 



We thank A. Bartl and W. Majerotto for many useful discussions. We are grateful to V. Latussek for his support in 
the development of the numerical program. This work was also supported by the Deutsche Forschungsgemeinschaft 
under contract Pr 1064/2-2 and the 'Fond zur Forderung der wissenschaftlichen Forschung' of Austria, Project No. 
P10843-PHY. 



[1] M.E. Peskin, Prog. Theor. Phys. Suppl. 123, 507 (1996). 
[2] J.L. Feng, M.J. Strasslcr, Phys. Rev. D 55. 1326 (1997). 

[3] A. Bartl, H. Fraas, W. Majerotto, B. Mosslacher, Z. Phys. C 55, 257 (1992). 
[4] S. Ambrosanio, B. Mele, Phys. Rev. D 52, 3900 (1995); 

S. Ambrosanio, B. Mele, Phys. Rev. D 53, 2541 (1996). 
[5] S.T. Petcov, Phys. Lett B139, 421 (1984). 

[6] S.M. Bilenky E.G. Christova, N.P. Ncdelcheva, Bulg. Jour, of Phys. 13,4 , 283 (1986). 
[7] Y. Kizukuri, N. Oshimo, Phys. Lett. B 249, 449 (1990). 

[8] S. Kawasaki, T. Shirafuji, S.Y. Tsai, Progress of Th. Ph., 49/5, 1656 (1973). 
[9] D.A. Dicus, E.C.G. Sudarshan, X. Tata, Phys. Lett. B 154 , 79 (1985). 
[10] S. Gottlieb, T.J. Weiler, Phys.Rev. D 32, 1119 (1985); 

S. Matsumoto, K. Tominaga, O. Terazawa, M. Biyajima, Int.J.Mod.Phys. A 8, 1153 (1993). 
[11] C. Dionisi, K. Fujii, S. Giagu, T. Tsukamoto, in Physics at LEP2, CERN 96-01, eds. G. AltareUi, T. Sjostrand, F. Zwirner, 
337. 

[12] H.E. Haber, G.L. Kane, Phys. Rep. 117, 75 (1985). 

[13] A. Bartl, H. Fraas, W. Majerotto, Phys. Rev. D40, 1594 (1989). 

[14] A. Bartl, H. Fraas, W. Majerotto, Nucl.Phys. B 278, 1 (1986). 

[15] H.E. Haber, in Proceedings of the 21st SLAG Summer Institute on Particle Physics, Stanford 1993. 
[16] H. Pilkuhn, The Interactions of Hadrons, North-Holland Publishing Company, Amsterdam 1967. 
[17] L.J. Hall, J. Polchinski, Phys. Lett. B 152, 335 (1985); 

A. Bartl, M. Dittmar, W. Majerotto, in Proceedings of the Workshop e+e" Collisions at 500 GeV, ed. P. Zerwas, DESY 

92-123B, 603 (1992). 

[18] A. Bartl et al., in Proceedings of the 1996 DPF/DPB Summer Study on New Directions for High-Energy Physics (Snowmass 
96). 

[19] Particle Data Group, Phys. Rev. D 54, 65 (1996). 
[20] S. Hesselbach, H. Fraas, Phys. Rev. D 55, 1343 (1997) 
[21] S. Ambrosanio, B. Mele, Phys. Rev. D 55, 1399 (1997). 



11 



FIG. 1. Feynman graphs for production and leptonic decay. 



FIG. 2. The amphtude squared for the combined process including spin correlations is composed by the spin density matrix 
for the production and the decay matrix times the pseudopropagator A2 squared (compare Eq. (9)). 

FIG. 3. Configuration of momenta in the laboratory system. The lepton angular distribution refers to 0m-, the opening 
angle distribution to Q-| — 

FIG. 4. Lepton angular distribution for ^/s = 193 GeV in scenario A for mo = 90 GeV with spin correlations fully taken 
into account (upper solid) and for assumed factorization (upper dotted); for mo = 200 GeV with spin correlations (lower solid) 
and for assumed factorization (lower dotted). 

FIG. 5. Lepton angular distribution for y/s = 500 GeV in scenario A for mo = 90 GeV with spin correlations fully talcen 
into account (upper solid) and for assumed factorization (upper dotted); for mo = 200 GeV with spin correlations (lower solid) 
and for assumed factorization (lower dotted). 

FIG. 6. Lepton angular distribution for ^ — 193 GoV in scenario B for mo = 90 GeV with spin correlations fully taken 
into account (solid) and for assumed factorization (dotted). 



FIG. 7. Lepton angular distribution for ^/s = 193 GeV in scenario B for mo = 200 GeV with spin correlations fully taken 
into account (solid) and for assumed factorization (dotted). 

FIG. 8. Lepton angular distribution for ^/s — 500 GoV in scenario B for mo = 90 GeV with spin correlations fully taken 
into account (upper solid) and for assumed factorization (upper dotted); for mo = 200 GeV with spin correlations (lower solid) 
and for assumed factorization (lower dotted). 

FIG. 9. Lepton angular distribution for ^/s — 193 GeV in scenario C for mo = 90 GeV with spin correlations fully taken 
into account (solid) and for assumed factorization (dotted). 



FIG. 10. Lepton angular distribution for ^/s = 500 GeV in scenario C for mo = 90 GeV with spin correlations fully taken 
into account (solid) and for assumed factorization (dotted). 



FIG. 11. Opening angle distribution in scenario A for = 193 GeV and mo = 90 GeV (upper solid) and mo = 200 GeV 
(lower solid); for ^/s = 500 GeV and mo = 90 GeV (upper dotted) and mo = 200 GeV (lower dotted);. 



FIG. 12. Opening angle distribution in scenario B for y/s = 193 GeV and mo = 90 GeV (upper solid) and mo = 200 GeV 
(lower solid); for ^/s = 500 GeV and mo = 90 GeV (upper dotted) and mo = 200 GeV (lower dotted). 



FIG. 13. Opening angle distribution in scenario C for ^/s = 193 GeV and mo = 90 GeV (solid); for ^/s = 500 GeV and 
mo = 90 GeV (dotted). 



FIG. 14. Energy distribution in scenario B for ^/s = 193 GeV and mo = 90 GeV (upper solid) and mo = 200 GeV (lower 
solid); for ^/s = 500 GeV and mo = 90 GeV (upper dotted) and mo = 200 GeV (lower dotted). 
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+1 


-1 
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+1 


+1 


-1 


-1 




P2 <-> P3 



TABLE I. Substitution rules for Wab,cd, see Sec. 3. 





tan/3 


M 




riirui 




X, X, 






<t 




-p90 


T^200 

x" 
x^ 


A 


2 


84 


-250 


46 


97 


97 


123 


104 


217 


207 


35.5 


1.33 


B 


2 


112 


448 


51 


98 


-97 


139 


110 


226 


210 


15.0 


0.74 


C 


2 


215 


-83 


76 


-109 


97 


214 


114 


279 


228 


23.6 


23.7 



TABLE IL Parameters M, /u and mass eigenvalues in GeV, total width of X2 in keV. The 
superscripts denote the value of the scalar mass mo- 





X? 


xl 




\ z \ H^^\ m) 


( U Z \ H^^ \ H^, ) 


A 


(+.94 1 -.33 1 -.08 1 -.08) 


(-.35 1 -.89 1 -.16 1 -.23) 


B 


(t.79 1 -.60 1 +.11 1 +.07) 


(-.62 1 -.76 1 +.17 1 +.10) 


C 


(-.17 1 +.22 1 -.19 1 +.94) 


(-.05 1 +.29 1 -.92 1 -.26) 




( B 1 1+3 1 1 ) 


( B 1 1 U'i 1 H'^ ) 


A 


(+.98 1 +.16 1 -.08 1 -.07) 


(+.12 1 -.95 1 -.16 1 -.23) 


B 


(+.98 1 -.15 1 -.13 1 +.09) 


(-.18 1 -.96 1 -.18 1 -.10) 


C 


(-.25 1 +.11 1 -.19 1 +.96) 


(-.18 1 +.23 1 +.92 1 +.24) 



TABLE III. Neutralino eigenstates. 
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^/i = 193 GeV 


= 500 GeV 




mo = 90 GeV 


mo = 200 GeV 


mo = 90 GeV 


mo = 200 GeV 


A 


38.9 


11.1 


32.7 


21.6 


B 


10.0 


2.5 


10.9 


6.1 


C 


23.3 


23.4 


6.0 


6.1 



TABLE IV. Total cross sections for e ^ X1X2 a-nd subsequent leptonic decay, X2 ~>' Xi'^^^ 
for = 193 GeV and ^ = 500 GeV with mo = 90 GeV and mo = 200 GeV. 



Afb = [cr(cos e_ > 0) - fr(cos 6- < 0)]/[ct(cos 9- > 0) + a(cos 9- < 0)] /[%] 




= 193 GeV 


= 500 GeV 




mo = 90 GeV 


mo = 200 GeV 


mo = 90 GeV 


mo = 200 GeV 


A 


10.2 


11.6 


2.3 


3.0 


B 


-0.8 


5.9 


-0.3 


1.8 


C 


0.07 


0.11 


-0.01 


0.00 



TABLE V. Forward-Backward- Asymmetry Ap g ior ^/s = 193 GeV and ^/s = 500 GeV with 
mo = 90 GeV and mo = 200 GeV. 
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PRD, Moortgat-Pick/Fraas, FIG. 1 
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PRD, Moortgat-Pick/Fraas, FIG. 2 
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PRD, Moortgat-Pick/Fraas, FIG. 3 
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PRD, Moortgat-Pick/Fraas, FIG. 6 
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PRD, Moortgat-Pick/Fraas, FIG. 8 
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PRD, Moortgat-Pick/Fraas, FIG. 9 
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PRD, Moortgat-Pick/Fraas, FIG. 11 
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